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STRUCTURES OF DOMINANCE RELATIONS 


RoBeErtT L. Davis* 
UNIVERSITY OF MICHIGAN 


Dominance relations and their siruciures have been studied in the context of animal sociology 
(Rapoport, 1949a, b; Landau, 1951a, b), and occur in many other theoretical models of the 
social and biological sciences. When Rapoport and Landau wrote, there was no method known 
for determining the number of dominance structures which could be defined on a set of n ele- 
ments, for » > 4. The answer to this question can be important in certain further investiga- 
tions of the structural properties of dominance. Using a general method developed elsewhere 
(Davis, 1953), this paper derives a formula to answer the question for any m, and gives an 
application of the method to treat a case not previously analyzed. 


1. Definitions and examples. A series of papers in this journal (Rapo- 
port, 1949a, b; Landau, 1951a, b) called attention to the importance in 
social science and elsewhere of a kind of relation seldom previously recog- 
nized, vz., one which was asymmetric and “connected” [in Russell’s 
(1919) sense] but not necessarily transitive. These were aptly termed 
dominance relations. 

To be definite, let V = {1,...,m} be any set of elements. Then we 
can describe any dyadic relation on WN (to itself) in terms of a matrix 
A = (a;;), where a;; = 1 just when 7 stands in the given relation to 7, and 
otherwise @;; = 0. 

Definition A. The dyadic relation A is a dominance relation if, and only 
if, (i) all a;; = 0, and (ii) for all 7 #7, either a;; = 1 or else a; = 1, but 
not both. 

Rapoport and Landau examined these relations in a model for the 
“neck-right” among barnyard fowl, also mentioning their applicability to 
other situations commonly obtaining in small groups. Any tournament, 
for instance, in which each participant meets another just once (no 
draws) gives rise to a dominance relation; this example drives home the 
point that the relation may, or may not, be transitive. The “Method of 
Paired Comparisons” in psychological scaling gives rise (in the modern 
view, where intransitivities are not simply thrown out as “‘error”’) precise- 
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ly to the structure of some dominance relation. And in genetics, the as- 
sumption of “dominance” in the one-locus diploid case amounts to assert- 
ing that there is some dominance relation defined on the set of alleles. 
Under the usual Hardy-Weinberg hypotheses, then, the relations between 
gene frequencies and phenotype frequencies are determined by the domi- 
nance structure. 

2. The structure problem. In such questions an investigator is seldom 
concerned with the particular assignment of names to the objects or indi- 
viduals being examined, but rather with the nature of some underlying 
structure. By structure is meant just what is left of the relation after the 
names are forgotten: the “map” (Russell, Joc. cit.), or graph, of the rela- 
tion. For a precise definition of structure, it is simplest to start with what 
it means for two relations to have the same structure. 

We want to say that two relations A and B defined on the set V have 
the same structure, or are isomorphic, whenever there is some permuta- 
tion + of N such that A has the same matrix with respect to NV as has B 
with respect to 7(NV). Technically, it is convenient to define, correspond- 
ing to each permutation 7 of the symmetric group S, on letters, a trans- 
formation ¢, which operates on our relations 


ty (A) - (Gx(a)xl3)) » where A = (a;;) ‘ (1) 


(Note that if P, is the permutation matrix corresponding to 7, then 
t,(A) can be thought of as P,A P=".) 

Definition B. Two relations A and B defined on N are isomorphic if, and 
only if, there is a permutation 7 in S, such that ¢,(A) = B. 

Definition C. The structure of the relation A is the class of all relations 
with which it is isomorphic. 

‘This is what Landau (1951a) called the dominance structure; he inserted 
the qualifier because Rapoport had used the word “structure” differently. 
To obtain a manageable set of invariants the latter (Rapoport, 1949a) had 
defined as “‘structure” something we may call the structure sequence of a re- 
lation A: this is just the set of row-sums of the entries of A. (These are 
usually written down asa sequence, in non-increasing order; this amounts 
to disregarding their original order.) From its definition, it is clear that 
this structure sequence is invariant under the group of transformations 
t,, 1.e., that isomorphic relations will always have the same structure se- 
quence. That the various such sequences do not form a “complete set’’ of 
invariants—in other words, that different structures may have the same 
structure sequence—can be seen already when m = 5 (Landau, 1951a). 

Now it is an immediate consequence of the definition that the number 
of dominance relations on a set of n elements is 2"-/2, But how many of 
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these are “essentially different” relations? That is, how many dominance 
structures can be defined on NW? 

The answer to this question may be useful in two ways: (1) The formu- 
la, though practically impossible to evaluate for sizeable n, may yet pro- 
vide a guide in, e.g., comparing proposed approximations or estimating the 
error consequent upon various assumptions; (2) Perhaps more important 
is that methods used here furnish tools for analyzing structural properties 
of dominance which could be applied, without really undue difficulty, to 
give a complete treatment of the cases, say,  < 12. The cases previously 
so treated—z < 4—were exceptional, while some of the more variegated 
aspects of dominance set in for m = 5 and n = 6. Finally, it is not unlikely 
that in certain theories where the number of “individuals” involved is 
small by hypothesis (as in some concerning social groups) these accessible 
cases m < 12 would suffice for useful applications. 

3. Algebraic background. Derivation of this formula rests on the theory 
of permutation groups and on a method developed (Davis, loc. cit.) to an- 
swer the more general question: How many non-isomorphic m-adic rela- 
tions (of all kinds) are defined on a set of m elements? It would be bootless 
here to repeat the development of that paper. What can, perhaps, be 
profitably undertaken is: (1) to state those results in sufficient detail to 
show their application to the dominance question; (2) on that basis, to 
derive the main formula in such a way that anyone familiar with the other 
paper may verify it; (3) to illustrate the methods in an example; and 
(4) to give the more readily obtained numerical answers. 

Suppose, then, that D, is the set of all dominance relations on NV. The 
symmetric group 5S, “‘acts on’’ D, in the sense that any transformation f,. 
corresponding to an element of S,, maps members of D, into other mem- 
bers of D,. Now, by an orbit in D, we will mean a set consisting of some 
dominance relation, A, and all its images, ¢,(A), as 7 ranges over Sy. 
Hence to count all structures in D, we have simply to count the orbits 
under S,,, since these orbits are precisely what Definition C says structures 
are. Further analysis now rests on the theorem that in any such case the 
number of orbits is equal to the average number (loosely speaking) of re- 
lations fixed per group element. More precisely, let s#(m) denote the num- 
ber of dominance structures on 7 elements and for each 7 in S,, let f(x) 
be the number of dominance relations A such that f,(A) = A. Then it can 
be shown that 


si(n) => i (n) (2) 


“TES, 


We can further observe that if + and ¢ are conjugate elements of S,, 
(ie., 6 = x— x for some x), then f(¢) = f (ar); this ‘says f is what is known 
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as a class function. Thus it is enough, for our purposes, to evaluate f on 
one representative of each conjugate class, [7], and then multiply by the 
number, c(7), of group elements in that class. The formula then becomes 


si(n) = D> 6 (x) f (n) (3) 
* [] 


where the summation is now over all conjugate classes, [7], in S,. 

In the symmetric group each conjugate class is determined by an - 
tuple of non-negative integers (f1, .. . , pn)—where fp; is the number of 
cycles of length % in the disjoint-cycle representation of any member of 
the class (consequently, 1: + 2p. +...-+ np, = m; the px’s define a 
partition of the integer m). And now, rather than the coefficient c(z), we 
will prefer to take the constant term inside the sum and thus use 


1 1 
n(x) = mi (r) = 1?1p,!2?2po!...m?np,! - 


[making use of the formula for c(z)], so that the latest version of the for- 
mula we seek is 


si(n) = >) n(x) f(x). (5) 


[x] 


Since for any » there is an explicit, though arduous, means of writing 
down all the partitions which determine our summands, and since we have 
determined n(7) for each of these summands in (4), it remains only to find 
an expression for f(a). This can be done by a method which is almost the 
same as those of theorems 4 and 5 (Davis, loc. cit.), deriving the number of 
structures for symmetric and asymmetric relations. The sketch which fol- 
lows is to be taken as a proof only in conjunction with the other paper; 
readers more interested in application of the formula might well skip to 
the final section, pausing only to observe the results of section 4. 

4. Derivation of the formula. There is one new wrinkle in the dominance 
case which makes later computations much easier. 

Lemma. If the disjoint-cycle representation of 7 contains any cycle of 
even length, then f(r) = 0. 

Proof. There is no loss of generality in assuming the given cycle to be 
[123 .. . (2k)]. Now by formula (1), #,(A) = A will require, among other 
things, that 


a a 


tht. 99 b¢a ++ = & op = 4) =: - 


But one of the entries a,44; or @x41,, must be zero and the other must be 
one for A to be a dominance relation, so they could not thus be equated. 
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We now may as well banish these even-cycle cases altogether, replacing 
the coefficient (7) defined in (4) by 


ai jae if + has any cycle of even length 
n (a) , otherwise . 


(6) 


As in several theorems of the former paper, the approach to evaluating 
f(x) is through a concept that may roughly be described as “the number 
of degrees of freedom in a matrix scheme if the matrix is to be fixed under 
i,.” That is, for each permutation 7 we define a number d(r) with the 
property that if you want to write down an arbitrary matrix with #,(A) 
= A, d(m) is the number of places in the matrix at which you can freely 
choose to write either a zero or a one. From this description of d(z), it is 
evident that our f(7) = 2%”. This gives us the final reduction 


si(n) = >> s(x) 24), (7) 


[x] 


Now it only remains to evaluate d(r). 
Theorem. The number of dominance structures on a set of x elements is 
given by formula (7) with 


d(r) = Ep -— 1 + Dabs, i) (8) 
i=1 


i<j 


where (7, 7) is the greatest common divisor. : 

Proof. The argument here—we are only dealing with 7 whose cycles are 
all of odd length—parallels those of the theorems mentioned above. Its 
formulation rests on our choice, in each conjugate class, of a representative 
a for which it will be specially easy to compute d(m), and further on our 
then splitting the matrix into “blocks” whose rows and columns will be 
permuted among themselves by individual cycles of z. 

The right-hand sum in (8), then, is precisely the same as that part of 
the corresponding formulas in the theorems about symmetric and asym- 
metric relations: in all three cases this part arises from those blocks whose 
rows and columns are acted on by cycles of different lengths. 

As for the left side: there are, for each i from 1 to n, p,(p; — 1) blocks 
(previously called “near-diagonal”) whose rows and columns are gov- 
erned by different cycles of the same length. In each of these, when we 
write out the requirement that #,(A) = A, we get 7 strings of equated 
entries (i entries per string), all of which can here be freely chosen so that 
A remains a dominance relation. But having made this set of choices for 
one such block, we have completely determined the choices in the “trans- 


136 ROBERT L. DAVIS 


pose block” (the reflection across the diagonal). Thus, from such blocks 
we get altogether (¢p;/2)(p: — 1) degrees of freedom. 

Finally, the “diagonal blocks” are those whose rows and columns are 
acted on by the same cycle. Here au, for instance, is zero, while of the re- 
maining entries in the first row we can freely choose just half before we 
completely determine the whole block. For i,(A) = A again gives 7 strings 
of eyuated entries, one consisting of diagonal terms and so permitting no 
choice, while of the remaining strings choice for any one determines that 
of its “transpose string.” Thus for these blocks we have p,(¢ — 1)/2 de- 
grees of freedom; the sum of these expressions for near-diagonal and diago- 
nal blocks then gives the left side of (8). A briefer and clearer form of (8) 
is easily seen to be 


dim) =) DS pbs D =D Ht. (9) 


Ce ee! 


5. An example. Actual counting can be systematized in a tabular set- 
ting. Since conjugate classes in S;, correspond one-one to partitions of the 
integer ”, we may list all these partitions in the first column [using the 
standard notation: e.g., (31?) stands for the partition of 5 into three parts 
given by 5 =.3 + 1-+ 1]. The second column repeats this information in 
terms of the non-zero cycle numbers. We next compute the denominator 
of the coefficient s(z) as given by formula (4), while the fourth column 
gives s(z) itself. The last column gives d(x); this is computed on the basis 
of (8) in the lines beneath Table I. Below these we then write down the 
sum giving st(m), as read off from formula (7) and this table. Note that 
the lemma above says we can discard all those partitions containing any 
even parts. 

To compromise between complete transparency and sheer tedium, we 
cite in illustration the case 2 = 5; this is the first which was not previous- 
ly analyzed by other means. 

Granted there are 12 dominance structures on a set of five elements, can 
we say what they are like? Now while the methods above provide the only 
way known to the author of finding sé(m), they are not well adapted to 
determining which relations lie in what orbits. However, the theorem ex- 
pressed in formula (2) rests on two hitherto unmentioned facts which 
simplify the task: (i) the set of all permutations leaving any relation, A, 
fixed is a group; (ii) the number of relations in the orbit of A is then given 
by the index of this subgroup in S,. 

Thus when » = 5 there are just three kinds of jix-groups, as we may 
call them: (i) those generated by 5-cycles, (ii) those generated by 3-cycles, 
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and (iii) the group consisting of the identity permutation alone. Here 
then, all orbits are of size 24, 40, or 120. (Note that members of the wes 
orbit all have conjugate—but not identical—fix-groups, and also that the 
same group may fix relations in different orbits.) 

The quick way now to a complete analysis of the structures is to use 
every available tool. In particular, the structure sequences described in sec- 
tion 2 afford valuable guideposts. Whenever we know that there is only 


TABLE I 

Partition Cycle-numbers ae Pel s(m) d(r) 
(S) ps=1 (S')1! 5 2 
(41) * * 0 * 
(32) x ~ 0 . 
(312) re eee (31)11(12)2! 1/6 4 
(221) * * 0 * 
Pei owes iS) | ot WS 
(15) jimi (15)5! 1/120 10 


[The sum is 
d(5) =$(5—1) +0=2, 


g(34%) =3t3—1) +3(2—1) +2=4, 
d (18) =§(5—1) +0=10, 
st(5) =2(4) +2(16) +73, (1024) =12.] 


one orbit with a certain structure sequence, we will know all about the 
orbit: just write down a relation (or graph) with that sequence to repre- 
sent the orbit and describe its structure. Now for 7 = 5 there are only 
nine structure sequences as against twelve orbits. However, combining 
structure-sequence arguments with those on fix-groups and orbits, we can 
derive Table IT. 

Thus the structure sequences are sufficient to determine seven of our 
twelve orbits unambiguously, and our description will be complete upon 
examination of the remaining cases. (Of course, this is the reverse of the 
investigating process; there you must first determine which are the am- 


biguous cases.) 
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Now, let 
O ih al ak Oa TiO 
() @ ab ah i ORO aera! 
A,=[0001 0] and As=]0001 1 
OF0R0RORI OT ORO eOmt 
fe Om te ORO 19.0.0: O20 


Inspection shows both these matrices have structure sequence (3, 3, 2, 
1, 1); further, since neither is fixed under any nontrivial permutation, 
both lie in 120-element orbits. That these must be different orbits is seen 
from the fact that a “man” dominating only one other dominates a “lead- 
er” in Ag, while this is not so for 47. Hence we can take A; and As to rep- 
resent the seventh and eighth orbits. 


If we further define 
OPFOR tae tet O+ Lod Ord O- Oat 20 
10100 ie tiiege ae Ue LOS OLE GO 
Aga=}00011 7], An=]? 000114, An=]} 0101047, 
ORT Oem T= O202 Uae 0.0" Ot 
Ot 0050 = FOS O80 dade OS O 


we are faced with three relations whose sequence is (3, 2, 2, 2, 1). It can 
be argued as above that A» and Ajo belong to 120-element orbits; An, on 


TABLE II 
came ar ta Re 
1 120 (4,3, 2, 1,0) 
We er ve rere 
3 40 | (4,2,2,2,0) | 
4 120 (45922: Ps 
5 40 (3,909.00) 
6 $2004. NS, ayasaoOen 
7 120 
A ian ($)352,.251) 
9 120 
ma hh gata 120 (3,2,.2.2.1) 
11 40 
12 24 (25.252) 292) 


ST a ey 
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the other hand, is fixed under the group {I, (123), (132)} and so belongs 
to a 40-element orbit. (Thus Ay; is more symmetric than the others in that 
all the “middle men” in Ay play similar roles.) That Ag and Ay) are further 
not isomorphic to each other is clearest, again, in simple combinatorial 
terms. In both the leader is dominated by one middle man; but in A that 
middle man picks up his other domination at the expense of a second 
middle man, whereas in Ajo the middle man dominating the leader also 
dominates the subordinate, and does not dominate either of the other 
middle men. These three relations, then, can be taken to represent the 
orbits 9, 10, and 11. 

Finally, Table III compares the number, Dom(m), of dominance rela- 


TABLE III 

n Dom(2) st(n) 

2 2 1 
3 8 Z 
4 64 4 
5 1,024 12 
6 32,768 56 
7 2,097,152 456 
8 268 ,435 ,456 6,880 


tions with st(m) for n < 8. These were quite easy pencil-and-paper com- 
putations, and despite the rapid increase of the number of partitions I be- 
lieve their number could easily be doubled with a desk computer. 
Landau’s inequality (1951a), 
Qn(n—1)/2 
Phen’ 


st(n) 2 


gives not only a lower bound but a very good approximation as soon as” 
n is moderately large. For example, for n = 8, st(8) = 6860, while 278! 
is 6658 to the nearest integer. From the viewpoint of this paper, this 
approximation means using only the partition (1"). If to this is added the 
term for the partition (31"-%), we have a much better approximation, 


namel 
ys Q(n?—n)/2 2 (n?—sn+ 8) /2 


n!\ (n— 313) ° 
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AN EXTENSION OF THE THEORY OF RESONANCES 
OF BIOLOGICAL CELLS: I. EFFECTS OF 
VISCOSITY AND COMPRESSIBILITY 


EUGENE ACKERMAN 
Puysics DEPARTMENT 
PENNSYLVANIA STATE UNIVERSITY 


The effect of the compressibility and the viscosity of the medium on the resonance of bio- 
logical cells is discussed mathematically. The resonant frequency is only slightly altered, but 
the motion is highly damped. 


The mathematical theory of mechanical resonances of biological cells 
was discussed in an earlier article in this journal (Ackerman, 1951). Two 
different cell models were considered, one a cell surrounded by a mem- 
brane with an interfacial tension, and the other a cell surrounded by a 
somewhat rigid cell wall. It was shown that either led to reasonable values 
for the physical constants of the cell membrane (or wall) computed from 
the observed resonances. However, the theory included several approxima- 
tions, which might limit its validity when applied to an actual cell. It also 
failed to explain why the resonances had not been observed by other 
methods and was unable to predict the relation between the height of the 
resonant peak and the amplitude of vibration. 

This is the first in a series of four papers designed to bring the mathe- 
matical theory into closer correspondence with the biophysical experi- 
- ments. In this paper, the effects of the viscosity and the compressibility 
of the media will be considered. Successive papers will discuss the cross- 
section presented to plane waves, the relation of the observed curves to 
the computed Q of the resonance, and, finally, the changes produced by 
the nonspherical shape of all real biological cells. (The Q of a resonant 
system is defined as the number of cycles necessary for the amplitude to 
be reduced by a factor of e*.) These discussions will all be limited to the 
interfacial tension model since it is the simplest to treat in detail. More 
complex analyses could be developed for the rigid shell model but seem 


somewhat fruitless at the present time. 
The effects of viscosity and compressibility are determined by the con- 
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stants of the cells and the surrounding medium (usually water). The coeffi- 
cient of viscosity, p:, of the interior of large cells has been estimated to be 
several times that of water, ie., 2 — 15 centipoises (Heilbrunn, 1952). 
This is the static value; it represents the force necessary to move two 
parallel planes a centimeter apart, with a velocity difference of one centi- 
meter per second. It should not be confused with the absorption of com- 
pressional waves in a liquid, which is due primarily to other effects. To 
consider a reasonable case, a value for p; of 10 centipoises, i.e., 10—* poises, 
has been chosen for this paper. The coefficient of viscosity, p., outside of 
the cell is assumed to be close to that of water, i.e., about 10-? poises. 
(Throughout the remainder of these papers the subscript ; refers to inside 
of the cell and , to outside of the cell.) 

The compressibility, 1/B, enters the equations through the velocity of 
compressional waves in the two media, c; and c.. To account for the 
absorption of plane waves, both c; and c, may be considered complex. 
[These velocities, ¢; and ¢,, are more precisely defined by equation (3).] 
Values for tissues and protein solutions (Carstenson ef a/., 1953) show that 
the real part of both velocities is about 1.5 X 10° cm./sec. and that the 
imaginary part is much smaller in absolute value than the real part. 

The inclusion of viscosity and compressibility complicates the equation 
of motion and the boundary conditions. The complete equation of motion 
for a real fluid is (Page, 1935) 


pu = —Vptinv (V-v) +7Vv. (1) 


In this # is the pressure, p the density, and v is the vector velocity. The 
dot above the term indicates differentiation with respect to time, ¢. The 
components of v in spherical coordinates are denoted by 2,, v%, and v. A 
complete solution for this equation is given by 

=—VE+VxA } 


p=BYV-v= —BV*®. 


and 


(2) 


The velocity potentials, @ and A, must be solutions of the wave equations 


_B +$inw 
p 


= ey; e (3) 


and 
A= (eV; (a)?= te, (4) 
p 
The potential ® represents the compressional wave and A the torsional 


wave. The components of A in spherical coordinates are labeled A,, Ag, 
and A,. The added absorption of a plane wave could be included by adding 
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an additional imaginary term to c. Little data exist to indicate the correct- 
ness of c*, the shear wave velocity, for viscous liquid solutions. 

Since two media are considered, two sets of functions are needed for ® 
and A. Within the cells, 6;, A,;, and their first derivatives vanish at the 
origin. In the surrounding liquid, ®, and A, behave, as r goes to infinity, 
as (1/r) e**. Solutions expressed in spherical coordinates, r, 0, and y, suit- 
able to satisfying boundary conditions at the spherical surface, r = a, 
have the form 


®: = DC, jn (kir) Pi (cos 6) ee r<a 

©. = =Dy{ jn (Ror) +n (kor) }Pi (cos 6) &Me™* =r Da, Se 

Ay =Zbmjm(Rir) P, (cos 6) ee r<a 

Ayo=Zjm{ jm (her) + %m (Ror) }Pm (cos 6) eV e™ (6) 
= z} _ per lp cas" 8) ee r>a, 

Ao; = Een jn (ker) Pi (cos 0) &™ e™* r<a 

Ato =Zfa{ jn (hor) + tn (Ror) }Pa (cos 6) e'™ e™" (7) 
a2} Je pr ees 0) a en 126, 

and 

Ay; =Zamjm (ker) Pi, (cos 0) puues r<ga, 

Ago =ZBm{ jm (ker) +m (Ker) }P'™ (cos 4) ee (8) 
~ Z| Ee BP eipl ela tee! 1 rib, 


where k = w/c and k* = w/c*. The constants w, Cr, Dn, bm, dm, Cny fay 
am, and B» are all to be determined by the boundary conditions. The func- 
tions 7, and m, are the spherical Bessel and Neumann functions respec- 
tively (Morse, 1948). The function P}(cos @) is the /th associated Legendre 
polynomial of degree 7. The approximations for the functions A,o, Ago, 
and A,, are possible since | %5’r?| >> 1 for all permissible values of r. If, in 
addition, the vibrations do not alter the position of the center of mass, 
then one chooses and m in the above formulae to represent even and 
odd numbers respectively. For a particular mode it will be necessary to 
have / the same in all the series. 

In addition, the functions represented in equations (5) through (8) 
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must satisfy the following six boundary conditions at the surface of the 
sphere, r = a: 
1) v, continuous 


a8; 1 ie dane ab, 
i 0 Ap Se 
Sjiarsiig 06 a ae ny ar 
Z- (9) 
1 fe) - 0 ee] 
paces miei Bison BAGS 
7 sin alae om el a 
2) v» continuous 
1 06; PeigAe 1a ito, 
7 06 + rsind dy 7 gp (1 Au) of Or 
(10) 
1 ee ee | 
rsin@ dp 7 ope ; 
3) w SURE 
1 0; 1 0A,;; ¥ 10%; 
kik area @ Ai) veiled amit oF 
7 sin 6 ow THOM ely r sin @ Or 
(11) 
1a 1 OAgo 
+o 5, (7 A oo) Sa Pay a6 ’ 
4) T,, continuous 
IP our. Ve; [- 0 Vro ° (=) 
lace a oe ie 12 
nls ay ya *)|= if a0” or\r : ia 
5) T,y continuous 
1 opty ie [: O Vro 0 4) 
= — : 13 
ni; 90° E 08 cay r ; = 
6) difference in T,, ar! the surface restoring force 
a 0 pi 0 Vr0 
(- pt Bi v-vit mn 28)—(— p48 v-v. +0. 2) 
(14) 


ze Dg ae 9 2 Dri 0? re 1 
= (sa) [a (si 6 a tame 0 dy? es 

In the above, T;, represents the components of the symmetrical stress 
tensor in the liquid acting on the surface of the sphere in the 6-direction. 
The interfacial tension is represented by 7; all functions are assumed to 
depend on time as e“*'. Condition 6) which involves time derivatives is 
used to determine the frequency and the damping of the vibration. This 
leaves two additional equations needed to specify the particular mode to 
within one arbitrary constant. We first split equation (9) into two equa- 
tions so that the compressional and rotational components of v, are inde- 
pendently continuous. This gives the two equations at r = a 


dd; dB, 
or or’ (43) 
1 0 P 0Aop; 1 0 Aa 
r sin AR a a ees |- r sin ‘ey ee ee ee |: (ie 
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The first, (15), assures that these modes are similar to those discussed in 
the simpler approximation. 
An additional equation can be chosen as 


A pare at ae (ele) 


Thus a particular mode has been chosen which is relatively simple to 
handle analytically. Another choice might lead to higher values of Q, but 
those of this mode are sufficient to account for the experimental curves. 

Without solving the equations explicitly, we note that v, will contain a 
sum of terms, at least one of which will involve (sin @)-, unless either 
both A, and Ay vanish [which would make it impossible to satisfy equa- 
tions (10) and (11)], or 7 > 1. And, by similar reasoning, equation (13) 
will contain terms in (sin @)—! unless / > 2. Since both 7, and v, must re- 
main finite at all points, the requirement / > 2 is immediately imposed. 
This states that the lowest mode is nondegenerate, @ depending on 
P3(cos @) e”*”. This conclusion agrees with the observations on paramecia 
(Ackerman, 1952b), but suggests that frequencies reported for mam- 
malian red blood cells (Ackerman, 1952a) were not the lowest resonant 
mode. Thus any finite viscosity limits the number of possible resonant 
modes. 

An approximate solution will now be found for the lowest resonant mode. 
The factors ee will be understood and not written explicitly. Since 
|k2a?| as well as |k3a?| are much less than one the following approxima- 
tions are valid 

@; = C2 jo (kir) P3 (cos 0) ~Cr’ sin’ 6, ai 
@, = Do{ jo (Ror) tins (Ror) }P2 (cos 0) ~ Dr * sin’ 6. 


To satisfy exactly equations (15), (16), and (10) through (14), infinite 
series will be needed for the components of A as well as an infinite series 
for ®. In each of these series, the largest term will be the lowest term. For 
this approximation, all coefficients except the lowest (i.e., C, D, e2, fa, as, 
and 3) will be neglected; A, will be eliminated exactly. In the case of 0, 
and 2,/r which appear in equation (14), only the contribution of the rota- 
tional modes to P? will be included. A more exact analysis would be ex- 
tremely tedious and would not alter the nature of the results. The general 
procedure will be to express the amplitudes of all the terms in equation 
(14) in terms of the constant, C, in equation (18). 

The first step is to substitute equation (18) into equation (15) and solve 


for D, givin 
OF £7, BINInE BeL iat 
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This value of D is used wherever D occurs throughout the remainder of 
this paper. 

Next, f may be found in terms of e, and C by using equation (1). For 
this purpose it is convenient to expand sin @ in a series of terms P?. We 
note that 0/dy = 2i, and substitute equations (7), (17), and (18) into 
(10). This is solved for fp giving 


4 eM 1. 2Cai—e .} Rij, (Ria) sue LSbie Sp (19) 


In a like manner, equation (11) is used to solve for 83. In this equation the 
expansion of (d/d@ sin? 6) in terms of P%, is used. This time (7), (17), and 
(18) are substituted into equation (11) and solved for §3, giving 


B,= —iae*[0.4Ca+ a} Hj, (Ba) syle i@) t]. 20) 


To find a relation between C, é, and a3 (16) is now used. Substituting 
in (7) and (8) and then using (19) and (20) to eliminate f. and #; a rela- 
' tion is found between as and ep. Since | 1/k2a 
written pee es 


ay| j, (Ba) +i zi Ba) t_1.2e,i j, (Ra) +4 2 Hy (a) 
(21) 

_ Car 

~ ea 


Another relation between a; and e, may be found by operating on equa- 
tions (12) and (13). Since »;>> 7. while » is continuous, these equations 
may be rewritten 


1 0v,; 0 Ve; pls 
rPEy wa az (= =e 
1 maha 0 Uw oe 


rsin 0 dp or\r 


Substituting in equations (9), (10), and (11) gives a pair of equations from 
which A, can be eliminated. When equations (7), (8), and (17) are next 
substituted into this equation, expressing all terms as 1, sin 0, or cos? 6, 
we have 


2Ca? [sin 6 (cos’@ — 1)] = —2ei[ (1 — cos” 6) Rk? a jy (ka 
+10 jo (Ria) ] —5as[ (4 cos* @—5 cos’ 6+1) ki a’ j3 is (kta) (22) 


+ (10 — 60 cos” 6+56 cos* 6) 73 (ka)]. 
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Equation (22) again cannot be solved exactly and shows the need for an 
infinite series. But one can find a first-order approximation by expanding 
both sides as a series in P!.; in particular, the left-hand is 


2Ca? [18 P!—4 Pi , 
Therefore both sides of equation (22) are multiplied by sin 6 (cos? 6 
—1)d(cos @) and integrated from cos 6 = —1 to cos @ = +1. This gives 
3Ca” = — 2.8 erik? a? jy (Rea (+20 esi jo (R70) 
—1.4ask7 a? 73 (ka) + (asijs (7a). 


(23) 


Although equations (21) and (23) are not solvable by a simple unique 
form, some insight into the relative sizes of e:, a3, and C may be obtained 
by considering actual numerical values. Choosing P. caudatum values 
again, and keeping 7. = 10-' poises, 7; = 10-* poises, p = 1 gives 


2 3 ie, 
Ki ret eee x2X10-#= 4; ka = 26/4 
2 . . z im, 
Bs a? =i IO x4 x 10-4 = 40; ka = 6.3 67/4 
tear 
pe 8° 


Reference to tables shows that 72(z), 73(z), and their derivatives are num- 
bers of the order of one in the range z = 2. Thus except for isolated values 
of a, equations (21) and (23) can be solved only by values of e, and as of 
the order of Ca?/10 or smaller. Equation (21) shows that a3 ~ 1.2 ey. 
Although this does not represent an exact determination, it is sufficient to 
set an upper limit on the effect of A on the resonant frequency and Q. In 
particular, these terms enter equation (14) through v, and v,/r. Writing 
explicitly, 7, at r = a, we have 


ead ; 


08; ion. a 
74= — + | 55 (sin 6 Ayi) ay 


‘Or ' rsin @ 


The P? coefficient is given by the expression 


G 
— 2Ca-+ terms of the order of 10° 


Similarly one finds for the P; component of v,/ratr=a 


0 Up; fe) oe aly 2 (4) 
or = —2C+5.5a3 = r 6.5 est ar v ’ 
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ar = +8C+5.5a3 = a (2) — 6. Sci ~ (2). 


Therefore the damping term in equation (14) has the form 


G 
in O Uri 0 Vr0 = —2C{;+4n,} + terms of the order of A 


nym aa 


Thus the contribution of the rotational motion to either the damping 
or the resonant frequency must be a second-order effect. The vector po- 
tential, A, is necessary to satisfy the boundary conditions but is unneces- 
sary to compute the Q or the natural frequencies of the cell. 

To make the solution exact one could use entire series for A» and Ay. 
This would leave a small term in v, proportional to Pj which could be 
cancelled by adding an equal but opposite term to ®, which would need 
another set of series for A, and Ay, etc. But all of these would be minor 
corrections to small corrections and will not be pursued here. 

Equation (14) may now be rewritten as 
2 (ni — 4 no) ° 8T 
ee 


1.7? — 3 


=0. 


Solving for w, we have 
= 8 a= ni t4no mt4no_ 1 
Oe poeta ATS ae) +4 74 Fas, =o, (1+35), 


where w, represents the resonant frequency. 


Choosing a = 0.7 X 107 cm., w. = 7.5 X 10® sec! shows that as 
before 


aah: 
T == wpa? = 4 dynes/cm., 


but now also one can rere a Q for the system as 


* a Wo 
= = nm+4ne perk 
Even if ni Were as low as »o, Q could not be higher than 7. 

This seems like a very low value of Q if one looks at the observed break- 
down curves (Ackerman, 1952b). However, a later paper will show that 
the shape of these curves depends on the geometry of the sound fields and 
the distribution of small nuclei, as well as on the computed value of Q. 

The results of this discussion can be summarized as follows: 

1. The effects of compressibility and radiation are negligible compared 
to those of viscosity. 
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2. The viscosity limits modes to those depending on e'’, where / > 2. 

3. Rotational (shear) velocity potentials, A, are necessary to satisfy 
the boundary conditions, but do not have a large effect on the resonant 
frequency, w., or the Q of the system. 

4. Viscosity does not have a large effect on frequency w., but strongly 
damps the vibration. 


SYMBOLS USED 
Y = spherical coordinates 
n = coefficient of viscosity 
B = bulk modulus 
¢ = compressional wave velocity 
c* = torsional wave velocity 
subscript ; = inside cell 
subscript . = outside cell 
p = density 
v = particle velocity 
p = pressure 
® = scalar velocity potential 
A = vector velocity potential 
a = cell radius 
k = w/c 
k? = w/c* 
T,e, etc. = components of symmetrical stress tensor 
jn = spherical Bessel function 
MN, = spherical Neumann function 
P! = associated Legendre polynomial 
w = 2x frequency (1 + 7/2Q) 
@ = 2x resonant frequency (undamped) 
T = interfacial tension 
Q = number of cycles for the amplitude to be reduced by a factor of e 
subscript » = odd integer 
subscript , = even integer 
t=J/-1 
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NOTE ON THE INTERPRETATION OF TRACER 
EXPERIMENTS IN BIOLOGICAL SYSTEMS 


H. D. LANDAHL 
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It is shown that the results of C. W. Sheppard and A. S. Householder (1951) can be modi- 
fied to include open systems which may or may not be in the steady state. Some conditions 
are given under which simple chemical reactions can be included. 


It is the purpose of the present note to discuss the interpretation of 
tracer experiments in simple systems which are either open or closed and 
which may or may not be in a steady state. This problem has been consid- 
ered in considerable detail by Sheppard and Householder (loc. cit.) for 
closed systems in the steady state. We wish here only to cast the problem 
in somewhat different terms, to show how the above authors’ results can 
be extended to open systems not necessarily in the steady state, and to 
discuss conditions under which simple chemical reactions can be taken 
into account. In their papers only tagged molecules were considered as 
having been measured. Here we shall suppose that untagged molecules 
have been measured by chemical methods. Biological problems which one 
may attempt to treat in such a manner are those involving material ex- 
change between blood and tissues, for example. More recently, some as- 
pects of this problem have been treated in a different manner (Schmidt, 
1952). 

Consider a collection of » compartments with barriers which prevent 
free flow, but which may allow diffusion of some substance S. Let there be 
n distinguishable forms of this substance. Assume that the mixing within 
each compartment is rapid as compared with the diffusion across the 
barriers. Consider compartments 7 and k. Let A jx, = Ax; be the area as- 
sociated with the exchange between j and &. Let the permeability in going 
from k to 7 be denoted by hj, while that in going from j to k be hz;. The 
values of hj, and /,; are equal in simple diffusion, in the absence of solu- 
bility effects, etc. In many cases only the product of hj, and Aj, can be 
determined. It is the product that enters into the expressions except in 
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certain ratios where the area cancels. Assume further that the flow rates 
are proportional to concentrations at least for a moderately large range 
of values about the normal levels. For the present assume that S is not 
involved in any chemical reactions. 

The flow from & to j of species r, whose concentration in k is "Ci, is 
given by A jx4j, "Cz, while that in the opposite direction is Axjhe; "Cj. 
Thus the net flow between these regions is the difference between these 
quantities. Summing over all the compartments , the net increase in the 
amount "S; of species r per unit time in compartment 7 of effective volume 
V; is given by 


as 1 r, r, 
a ('C;V;) = DS Ahir Cy — A jeltaz C3] - (1) 
k 


If the values of *C; are measured and also the V,’s, then the values of A,; 
and h;, are found from (cf. Sheppard and Householder, Joc. cit.) 


NG er (IC; V;) pres 
eds Aut i Oe n a A n, , 5 

A jehjx = a nk 
mC, : / nC, 


In a closed system in which no ’s are zero, one has att = ~ 


hiny "Cx (@) _ Choo 
ee ES WS Oe (3) 


Also we have, if "A; is the absolute specific activity of species r or fraction 


of all types of species which are of type r (Sheppard and Householder, 
loc. cit.), 


C= "As DC =" ALC; . (4) 
Hence (1) may be written as 
o's. ‘ 
“at 7 => Pe [ A; (hin AnCe ) a "A; (hiss pe c?)] : (5) 
jo 


We may write (5) as 
aS; = (TALS jx = "A ;dS4;) ) (6) 
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[which is equation (4) of Sheppard and Householder] with 


DS ee od Sah An Ce 
Sn Cae 15u = hn AnCe dl. (7) 

Now if each of the values C? remains nearly constant for a sufficiently 
long time, then the values dS; may be treated as constants and (6) may 
be solved. If each "C; is very small compared with "C;, e.g., all tracers are 
at low enough concentrations, then C? ~ "C; and hence in this case the 
system (6) can be solved if the values of "C; for various & are sufficiently 
constant. If they are in a steady state, then CZ/C? = Cyo/Cjo and thus 
DS jx = AS4; = NjpAyjsCEdt = hy;A4jsC?dt. The solution (2) becomes equiv- 
alent to (6a) of Sheppard and Householder if each row of both numerator 
and denominator is divided by C?, since "A; = ’C;/C?. If the above 
conditions do not hold, then it is necessary to return to equation (2). 
Note, however, that each of the concentrations of the various species 
may be measured in arbitrary units. This is evident from an inspection 
of either equation (1) or (2). 

We have assumed that S is not involved in any chemical reactions, al- 
though the possibility that these may be taken into account has been sug- 
gested by Sheppard and Householder. However, if S is involved in a re- 
action of the form 


Ry 
X+5S2 XS (8) 


in, for example, the mth compartment, then, if ’C,41 is the concentration 
of rth species of XS, we have 


OC a 41 


dl = Te. Ore a k-1'Ch41 (9) 


Vn 


where, if X and XS cannot leave the compartment, X = Xtotal — Cra 


If for any reason Xtotai >> C241, then the m compartment case can be 
changed to an x + 1 compartment case if Vay = Vn, MX = Intijn, ka 
=hp, x11. If &X varies appreciably but slowly enough, one would expect to 
find the calculated values of 4X = /tny1,n to vary slowly, with the other 
calculated results being relatively unperturbed. Note that if all the spe- 
cies except one are at very low concentrations and if this one is approxi- 
mately in a steady state, and hence nearly constant, X is approximately 
constant and the reactions of the type (8) may be included as additional 
compartments. 

If there is a reaction similar to (8), but in which XS breaks down at 
rate k, to form P and X, and if the product P is counted as S in an addi- 
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tional compartment, this type of reaction can be included if the condi- 
tions given in the previous paragraph hold. Then Ans, ny2 = 0, and 
Iin42,n41 = ko, other k’s being zero. If the last step is reversible, the rate 
constant being k_2, then nasi, nyo = RX. 


The author is indebted to Drs. C. W. Sheppard, A. S. Householder, and 
G. W. Schmidt for comments and suggestions. 
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The probability that one molecule will collide with another is calculated for the case in 
which the molecule is confined to equally spaced rings about the second larger molecule except 
when it jumps from one ring to the next. The result is practically the same as that obtained 
for an infinitesimal mean free path, not having first-order terms in the ratio of the mean free 
path to the particle radius as does the result of the calculations of R. Wijsman (1952). 


In a number of problems of biological interest a knowledge of the prob- 
ability of two molecules coming together is desirable. This problem has 
been considered by Wijsman (Joc. cit.). It is the purpose here to treat the 
problem in a different way and for a somewhat different case. For con- 
venience we shall refer to the larger molecule as a particle and suppose that 
it does not move. To calculate the probability of collision when the dis- 
tance between the particle and molecule is 71, we suppose that there are 
many molecules produced at the surface of a sphere of radius r = », and 
then calculate the fraction which reaches the particle at the center (Wijs- 
man, Joc. cit.). 

Consider a particle of radius r} which is a perfect absorber located at 
the center of a sphere of radius 7, at whose surface o molecules are pro- 
duced per cm.? per sec. We wish to determine the fraction of molecules 
which contact the particle so that the center of the molecule is a distance 
ro equal to rj plus the molecule’s radius. Let the molecules have a mean 
free path for scattering (/), and let the probability that a molecule be 
absorbed in unit time be a. We consider the molecules of the solvent 
(water) to be in concentric shells about 7» such that the centers are at 
ro + (wn + 4)1, being an integer. Hence the molecules may be thought to 
squeeze between the molecules of the solvent from 7 + (n + 1)1 to 
ro + nl so that the mean free path of the molecules is of the order of the 
diameter of the solvent molecules. To simplify matters we take / to be 


equal to this diameter. 
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Consider a molecule at 7) + ml. Let it be constrained to remain a dis- 
tance ry + ml + I after one step. Let P; be the probability that a mole- 
cule on the right of the particle moves to the left toward the particle, and 
let P, be the probability that it moves to the right. We shall suppose that 
the molecules move a distance / in one of six ways: 1) to the left to ro 
+ nl — 1,2) to the right to 7) + 1 + J, 3) to 6) up, down, in, or out but such 
that the distance is 7) + ml. Since the last four motions, which are equally 
probable, are not in the same plane, these would result in a net transport 
unless P, ~ P;. Hence let P, = 1/6+ «, P; = 1/6 — ¢, so that P, 
+ P; = 1/3. Determine ¢ such that the center of gravity after motion is, on 
the average, the same as before. If 6 is the projection toward the left of 
any one of the four perpendicular motions, each motion being for a dis- 
tance J, then one finds e = 6/3/. But from the geometry of the situation 
6/1 is equal to 1/2(ro + nl) if 1 < ro + nl, which always holds if we assume 


l < ro. Define 


eee eT eee (1) 


To 
Hence we have for the probability of moving to the right when starting 
in ring 2 
fa ee ee 


t93b ne Obert bee (2) 


Similarly P; = 1/6 — € = dn1/6dn. 

Let V,, be the number of molecules at 7) + nl. Let A, be the area of the 
nth ring so that A, = mrid;z. Let r be an average time between effective 
collisons. Then the number of molecules which move from ” — 1 to » in 
the time 7 is equal to the product of the number of molecules in the 
— lst ring and the probability that they move to the right, i.e., 


iy Gn 
Na-i e 1 Na-1 6¢n—1 . 


Similarly the number moving from + 1 to m is given by N41 times 


Piuhiy NatiPu nti= Nazi sc ; 
n+1 


the number moving away from ” to nm + 1 and to  — 1 is likewise found 
to be 


diye Packt pd nos Mint 

Sisnn ait acts dig 
Subtracting this last quantity from the sum of the former two quantities 
leaves an amount which divided by + is the rate of disappearance of mole- 


A 
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cules from ring 7. The rate of disappearance must just be that due to the 
reaction and equal to aN. Thus we find, if 


oT, (3) 
Yn+1 — (2 a7 6B?) SS ea es = 0 5) (4) 

where 
(es 6ra= or, (5) 


the quantity D being the diffusion coefficient given by /?/6r. 
If we introduce m = (7; — ro)/l and 


2 4 
w= 14+ Fy yieng (6) 
then the solution of (4) is 


(7) 


y,= Aw +Bo", n<m, 
yn = Fo" +Gu ”, n>™m. 


Clearly, since w > 1, F = 0. Now at m, where r = , N%, = Ni,, and 
hence y%, = yi, while Vj; = yi = 0 (rY = 7%). Furthermore, the inflow at 
n = m must equal the amount flowing away from the boundary or 


A lh mo 


— nai t (248°) ¥n— I n-1 = = (8) 
Hence 
i : dnAo mal? rs es 
Nn = VnOn = 4 (w FS): ie 
2Da"V\ e+e 
Pd, ( (9) 
e e Avhmo ti, (Ge eye Ome 
Nn= VnPn = ri 0 
2DV6:+5 
Define P,, by the relation 
DN, 
aoe (10) 
Pm ¥ Ag? a’ 


ie., the ratio of the number of molecules absorbed on the particle per unit 
time to the total produced per unit time at 7. Then substituting (9) in 
(10) we find 

P= ———. Cia) 


158 H. D. LANDAHL 


If B <1, and since ¢n = ¢m/¢o = 1/70, and m = (r, — ro)/l, then the 
probability P(r) that a molecule produced at a distance 7, from a given 
particle will collide with the particle is given by 


Pm=P(ri) = * e—mB (1—(1/24)8' +...) ~ Yo e—V4/D(r,—15) . (12) 
1 


ry 
This last expression does not contain a correction term in //ro as does the 
corresponding expression of Wijsman (Joc. cit.), otherwise it is the same. 


The author is indebted to Drs. G. Schmidt and C. S. Patlak for reading 
and discussing the manuscript. 

This investigation was supported by a research grant G-3312 from The 
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Health Service. 

LITERATURE 


Wijsman, R. A. 1952. “Probability that Two Diffusing Molecules Will Collide.” Bull. Math. 
Biophysics, 14, 121-26. 
RECEIVED 6-19-53 


—_—s?* 


eS 


BULLETIN OF 
MATHEMATICAL BIOPHYSICS 
VOLUME 16, 1954 


THE RESONANCE OF THE ARTERIAL SYSTEM 


GEORGE KARREMAN 
COMMITTEE ON MATHEMATICAL BIOLOGY 
THE UNIVERSITY OF CHICAGO 


The arterial system is assumed to consist of two elastic chambers connected by a conduct- 
ing channel. It is assumed that a current of fluid enters one chamber, whereas the other cham- 
ber is drained by a pipe with a certain peripheral resistance. The continuity of the fluid is de- 
scribed by a differential equation for each chamber. The inertia resistance of the conducting 
channel is taken into consideration. 

It is shown that the system may possess a resonance frequency. The latter, if it exists, as 
well as the damping coefficients are expressed in terms of the elastic moduli of the chambers, 
the conductivity of the channel, and the peripheral resistance. It is shown that with plausible 
values of the latter variables the resonance frequency as determined theoretically has the right 
order of magnitude as found experimentally. 


W. F. Hamilton and P. Dow (1938) discovered standing waves in the 
arterial system. These investigators (loc. cit., also Hamilton, 1944) as- 
cribed the production of these waves to the reflection of the advancing 
pressure wave in the arterioles and superposition of the refiected waves 
and of the incident wave. This explanation of the standing waves in the 
arterial sytem has been criticized recently on theoretical grounds by the 
author (Karreman, 1952) and on experimental evidence by R. S. Alexan- 
der (1953). The latter author (1952, Joc. cit.) has presented considerable 
evidence that the more or less sudden lowering of resistance at the site of 
the diaphragm where the major vessels, celiac, superior mesenteric, and 
renal arteries, branch off from the aorta plays a major role in the genesis 
of the aortic standing wave. The same author (1953) considers as the basis 
of the standing wave oscillation the fact that the retrograde transmission 
of the initial surge at the diaphragm back to the root of the aorta requires 
the same time as for the forward pulse wave to reach the resistance of the 
femoral bed. He (1953) therefore suggested considering the arterial system 
as consisting of two “closed” segments with equivalent wave lengths 
capable of producing a resonant oscillation about a diaphragmatic node, 
representing a to and fro surge of the aortic blood column between the 
root of the aorta and the femoral bed. 

We shall study in this paper a model of the arterial system which con- 
sists of two elastic chambers connected by a conducting channel. This will 
be an extension of the well known air chamber model of O. Frank (1899), 
in which the whole arterial system is taken to be an elastic chamber 
drained by the peripheral resistance of the arterioles. Frank himself (Joc. 
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cit.) also studied a two-chamber model, but one of the chambers was sup- 
posed to represent the venous system connected to the arterial chamber 
by a Poiseuille resistance representing the connecting arterioles and capil- 
laries. We shall suppose, as in the one-chamber theory, that the elastic 
properties of each chamber are the same everywhere in the chamber and 
that each chamber is so small that the pressure in it may be taken the 
same throughout the whole chamber. Let ©, G2, Pi, P2, denote the elastic 
moduli dP/dV and the pressures of the first and second chambers respec- 
tively (Figs.1A, B), 2 the current into the first chamber, 7 that from the 
first chamber into the second through the connecting channel BC, and w 
the Poiseuille resistance of the pipe #G draining the second chamber rep- 
resenting the peripheral resistance. The equation of continuity for the 
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first chamber is then given by 


1 dP, 


G dp ota (1) 


in which (1/G,)(dP,/di) represents the time rate of the amount of the 
elastically stored fluid, which together with the outgoing current i, bal- 
ances the incoming current 7. Similarly the equation of continuity of the 
fluid in the second chamber is represented by 


1 dP, , P»—P; 
©. di w 


=i, (2) 


expressing the balance of the incoming current i, with the storage term 
(1/G2)(dP2/dt) and the current (P; — P3)/(w) through the Poiseuille 
resistance w at the ends of which the pressures are P, and P3. Equation 
(2) is the equation of the one-chamber theory of D. M. Gomez (1937). 

To treat the fluid motion in the connecting channel BC we will proceed 
as follows. To illustrate the method we will assume here that the channel 
BC is cylindrical with radius r. The force on the fluid in the channel is then 
mr?(P; — P2), the mass of the fluid in the channel zr%p, so that Newton’s 
equation reads 


—- =—(P:—P,) , (3) 
p 


with c = (zr?/l), since diz/dt = zr?(dv)/(dt) in which dv/dt is the accelera- 
tion of the fluid in the channel. Equation (3) is similar to one of the fun- 
damental equations of the theory of multiple resonators as developed by 
Lord Rayleigh (1896, p. 173). The conductivity of the channel is c, the 
inverse of which is the inertia resistance of the fluid mass in the connect- 
ing channel. Ph. Broemser (1932) introduced such an inertia resistance in 
an extension of Frank’s studies of the air chamber type. However, he also 
studied only a one-chamber model with an inertia resistance in the en- 
trance lead. An extensive critical review of these and similar theories of 
the circulatory system is given by A. Apéria (1940). 
Substitution of 7 from (2) into (1) leads to 
ledPiyeloid Ps a Poe (4) 


aaa ae 


Elimination of P; from equations (3) and (4) gives 


p @P,, p_ @P2 (jt . Date Pa (5) 
cG,G, dé Giw de ae Gi As 
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To find the natural frequencies of the ie we have to solve 


p d?P, p d?P» (=+ dP» aly (6) 
cGi © dé G,cw d? = G! & Ap 
To achieve this we assume that 
Py = eit, (7) 


On substitution for P2 from equation (7) into equation (6) we obtain after 
some rearranging 


G» 


a+ — o? “(G + G:) a+———=0. (8) 


c “: G 

w 
Instead of solving equation (8) analytically, which would give rather un- 
wieldy expressions, we will first estimate the parameters p, G, G:, w, and 
c. The parameter p is of the order of 1 g. cm.-*. As introduced above, 
& = dP/dV. For a tube with radius r whose wall has an elastic modulus 


E and a thickness 6 we have (Karreman, Joc. cit., p. 335) 
2FT 
a4V =.V Es dP. 
Therefore 
G= dP E51 __ pa 
dV ..dr.¥ ay 


in which a is the velocity of propagation of pressure waves in such a tube 
given by the formula 


2 
Gat and Go. (9) 
As plausible average values we may take (R. S. Alexander, personal com- 
munication): a, = 500 cm. sec.—', a2 = 1000 cm. sec; i in Figure 1B, 
for which we will take the distance from the top of the arch to the dia- 
phragm, is about 18 cm. For V; and V2 we need also the values of r 
= 0.60 cm. and R = 0.72 cm., which values are derived from the figures for 
the cross-sectional areas of the lower thoracic aorta (1.15 cm.2) and of the 
abdominal aorta including the celiac, superior mesenteric, and renal 
arteries (1.62 cm.?) as given by Alexander (1953). Finally we have taken 
50 cm. as a plausible value for 2. This gives Vi = 20.3 cm.? and Vy 
= 81.6cm.* Substituting these values into (9) together with p = 1 g.cm 


Ne 
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gives 
©i= €.=1.23 X10! g. cm.—! sec.-2 . (10) 


Using the same values of a), a2, 7, and R, we obtain (E,/E») (51/82) = 0.3; 
therefore assuming that 5; ~ &, we have E; = 0.3, showing that the 
wall of the aorta between the arch and the diaphragm is more elastic than 
that below the diaphragm. To obtain the order of magnitude of w, we will 
suppose that all arterioles are parallel to each other. We have then for 
their total resistance w (assuming Poiseuille flow in them) 

Syl 1 


Ww — — 
TR* n’ 


in which 7 represents the viscosity of the blood, / the length of the ar- 
terioles, R their radius, and their number. Using the data compiled by 
H. D. Green (1950), we have R = 10 cm., » = 4.10’, n = 2.107 g. 
cm.— sec.—. Taking / = 1 cm., we find w = 10° g. cm.~* sec.—!, whereas 
1 = 0.1 cm. gives w = 10? g. cm. sec.—. 

To determine the order of magnitude of ¢ we remark that c = 1/o where 
o is the inertia resistance of the channel. To determine o we replace the 
channel by a series of connecting cylinders each of length dx and whose 
radii progressively increase from r to R. If 7 is the radius of such a cylinder 
the latter has an inertia resistance dx/a7”. Therefore 


d 
¢=2 = : 
aT. 
the summation being extended over all cylinders. As a crude approxima- 


tion we obtain in this way . 
saree 
st: focee (11) 
ae 


L being the length of the connecting channel. For 7 we will take the linear 

increase with x depicted in Figure 1B. 

1s eg 
L 


eairck- x. (12) 
Substitution of (12) for 7 into (11) leads to 


7 , 
ort (a3 
c=. (13) 


From what has been said above we then obtain 


cae. (14) 
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Substitution of the above-mentioned values of r and R, and taking the 
plausible value L = 5 cm., we see that c is of the order of magnitude of 
10-1 cm. It is hardly worth while to determine c more accurately due to 
the large uncertainties and our crude treatment, so we will keep only its 
order of magnitude. 

Using the given values of our parameters in (8) we see that one root a 
of that equation is negative and approximately given by 


so = 


= (G+ G) a+ = 0 
or 
2 
WwW 
eo ia eta (15) 
Ei,  & 


The two remaining roots a; and a3 approximately satisfy the equation 


1 
ab e/a Be pe “(G+ G) =0, (16) 
Wi deal atl 
& ©, 
whose solutions are 
- -3 gunsy ny 
2.3 = S Set iV oe) gets G.). (17) 


It turns out that with the values given above for the parameters © 
= & = 1.23 X 10* g. cm. sec, w = 10° g. cm. sec, c = 107 cm., 
and p = 1 g. cm.-%, the expression under the radical in (17) is negative 
and we have damped oscillations. However, a two- to three-fold increase 
of the radius of the arterioles would increase w 24 = 32-fold and would 
make the expression under the radical in (17) positive and therefore 
obliterate the oscillations. This might explain the disappearance of the 
standing wave system after the injection of a vasodilator as has been ob- 
served by Hamilton (1944). For an order of magnitude of the parameters 


as just mentioned, the frequency » of the oscillations is determined by 
the formula 


w= VE( G+ 2) . (18) 


Inserting the mentioned values of G, G:, c, and p we find 


wo = 50 sec.—!, (19) 


>” =. 


| 
: 
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From (19) we obtain » = wo/2m ~ 8 sec.-1, yielding for the period T of 
the oscillations 


i= 1-2) msec . (26) 


This compared rather favorably with the experimental finding of a T of 
the order of 140 msec. (Alexander, 1952). 

Further we will derive the solution of equation (5). Denoting, as before, 
the roots of equation (8) by a1, a2, and a3 and remembering that we have 
assumed P; constant in time, we can write (5) as follows 


eee eb-. ak) (Den) (2?) = (21) 


in which the symbol D represents the differential operator d/dt. To obtain 
the solution of (21) we have to add to the solution of the homogeneous 
equation corresponding to (21), 


(D+ ai) (D+ az) (D+ as) (P,—P;3) =0, (22) 


a particular integral of equation (21). The solution of (22) is given by 
P,—P3= Aeu!+Be* sin (wot + ¢) , (23) 


in which 8 is the real part Re(ae, a3) of ae or a3, wo the frequency given by 
(18), A, B, and ¢ arbitrary constants to be determined by the initial or 
final conditions. To illustrate the method we will assume that 


=I sinwt. (24) 


This is one of the two assumptions Frank (Joc. cit.) made about the time 
course of the current i; in his air chamber theory. Because 7, = 0 at the 
beginning and at the end of the systole we will assume that the systole 
starts at ¢ = 0, that it has a duration S, and that wS = 7 or 

(25) 


ao = 


YAIA 
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The particular integral of equation (21), after (24) has been substituted 
in it, is (Ince, 1944, p. 139) 


1 c&,C, 
Ep a ff t 
ie AAT OIE SA MONE CP UR gg 
c G2 1 9 
= a unniiD Asta el Datataal ea 
c&, E> (—D+ai1) (— D+ a2) (— D+ a3) : 
<p (= Dead) (= DF) (= Daz) 
2 
p 
ae DPS (ay+ao+ a3) VBS (ayao+aja3 + @2a3) D-+aja2a3 
s (Fat) (o®b a) (oP a2) 
Xsin wi 
(26) 
ee Spe £64 6.) D+ < = 
2 12 i 
<a Ate ath Gea cee 
- c€,C, 
p 
jo <(G.4 &,) @ Cos ot — S(t £51) sin wt 
a (a? a2) (a? Fal) (a? fal) 
6 
(w? — w?) w cos wt — —= (w? —w? ) sin wi 
_ ¢&,&, 27 
p (w? + at) (w?+ a2) (w+ a2) ; 
in which 
a a Es 27 
§=— (27) 


and w is given by (18). 

The solution of (21) is therefore the sum of the right-hand terms of 
(23) and (26) after, in the latter, the limits ¢ and 0 for ¢ have been substi- 
tuted and the corresponding difference has been taken. Transferring P3 
to the right-hand side we obtain for P» 


P.= A eu'+Be* sin (wot +) 


c GG, (w? — w5) @ (cos wl — 1) — S2(2 = sin wl (28) 
————] 


Tp | at bay Ga bay 
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During the diastole J = 0 and we have with other constants F, G, and 
y instead of A, B, and ¢ respectively 


P2 = Fev(*S) + G e(*-S) sin {wo (¢-— S) +h} +P. (29) 


To distinguish the P, given by (28) from that given by (29) we will denote 
the former by Ps and the latter by Ps». To determine the constants 
which have been left arbitrary so far we will assume that the diastolic 
pressure Pp, the pressure Ps at the end of the systole, and the slope of the 


curve (dP/dt),-s at that moment are fixed. Using (25), (28), and (29) we 
then have 


Pass =A+Bsin¢@+P;=Pp, (30) 
Pasyig= Ae -+Bo sin (pS +4) — SEAS 
2w (w” — wo) ie (31) 
* 1 (aA a2) (uo? Fal) (ot Fag) PsP a 
(Se) = Aare" + B&5{ B sin (woS +4) +9 cos (woS +4) } 
Se a ir (32) 
c&; &, pe or )e -($), 
Se TRE Ee ea 
Pop,-g =F +G siny +P; = Ps, (33) 
Papz. 4 =Fes12+Ge? sin(wD+y) +P;=Pp, (34) 


in which D = T — S is the duration of diastole and 


dP ; _ (dP 
ew) = Fa tG(Bsiny tuo cos y) <(F a) 


Because in all practical cases 8 = Re(a2, a3) is much smaller than », 
the factors w? + a? and w?+ aj in the denominators of the fractions in 
the right-hand members of (31) and (32) may be approximated by 
w? — w?. The fraction in the right-hand member of (31) contains, then, 
the factor 1/(w? — w?) which becomes infinite for w = wo. This shows 
explicitly the resonance property of the system considered here. 

The amplitude J is related to the stroke volume (SV) by the relation 


svy=f “I sin ot =F. (36) 


To illustrate the method we will suppose that the stroke volume is 65 
cm.’, the duration of systole 250 msec., that of diastole 600 msec., so that 
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the total duration of the cardiac cycle is 850 msec. and therefore the fre- 
quency of the heartbeat is 70.6 sec.—'. All these figures are about the aver- 
age values for the human being (Best and Taylor, 1945). With this value 
0.25 sec. for S we find w = 12.6 sec. from (25). For the length Z, in Fig- 
ure 1 we will take 36 cm., which is about the distance from the top of the 
aortic arch to the diaphragm. Assuming that a; = 500 cm./sec. we find 
288 msec. for the double round trip time to the diaphragm, which is about 
equal to the duration T of the standing wave (Alexander, Joc. cit.). There- 
fore T = 288 msec. and wy = 21.8 sec.—'. Assuming that the average cross- 
section of the thoracic aorta is 2.5 cm.? (Bazett et al., 1935) we find for 
the volume V; of the part of the aorta-above the diaphragm 36 X 2.5 
=90 cm.*. To find the volume V2 of the rest of the arterial system we will 
subtract the latter value from the total volume V of the arterial system. 
This procedure for the calculation of V2 is different from that in the be- 
ginning of the paper, but has the advantage that it does take into account 
the arteries which branch off from the major arteries and also contribute 
to the elastic modulus of the second arterial chamber. For the value of V 
we will use the figure 380 cm.’ which Apéria (Joc. cit., p. 181) calculates as 
the value of the volume of Frank’s air chamber for a particular human 
being, L.A; it should be mentioned here that he calculates in his own 
method a value two to three times larger. Though admittedly there is some 
uncertainty about this value we will keep for our illustration here the 
lower value. For the value of a2 we will take 750 cm./sec. as a reasonable 
value for the whole part of the arterial system below the diaphragm (cf. 
Bazett et al., loc. cit.). We find in this way according to our formula 


2 2 

@ =F with p=1 gramcm.°, ga oe 2780 dyne cm.-, 
2 

Beat Ce a 1940 dyne cm.~>. 


Substituting these values of ©, and ©, together with the above-mentioned 
value of wo into (18) we find ¢ = 0.1 cm., which is in agreement with our 
estimates of c above. Finally we will assume w = 100 dyne cm.~ sec. This 
value of w being probably at the low side, it might be remarked that it is 
possible that P is already constant after the flow has passed through a 
fraction of the Poiseuille resistance (Apéria, Joc. cit., p. 98). It is then 
necessary that Ps be sufficiently high. This is the case in (37) and (38). 
We find that after calculating the values of A, B, $, F, G, and wy with 
P; = 79 mm. Hg for the systolic part of the pressure P, in mm. Hg and 
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t in seconds 


Pg = 57 e—"-4#-4 57.7 e—4-0¢ sin (21.844 4.95) 


3 
— 22.9{ (cos 12.6¢—1) — 0.58 sin 12.6¢} +79 , sig 


and for the diastolic part of the pressure P» 

Py = 18.4 e—-4*- 14.0 e—4-% sin (21,842.01) +79. (38) 
The pressure P2 given as function of time ¢ for 0 < t < S by (37) and for 
S <i<T by (38) is represented in Figure 2. 
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A recently introduced approximation method is applied in order to obtain an expression for 
the amount of a substance remaining within a nerve at any time, the nerve having been soaked 
for a long time in a solution containing the substance until the time zero when it is bathed in 
the same solution but without the substance. The case of a uniform nerve without a sheath 
leads to substantially the same results as previously obtained by A. V. Hill (1928) for this case. 
A solution is given for the case of a nerve without sheath but having fibers which are per- 
meable. In this case it is shown how an effective diffusion coefficient for the interstitial fluid can 
be obtained, as well as the effective inward and outward fiber permeabilities. A solution is 
given for the case of a nerve with a sheath in which the substance considered does not penetrate 
the fibers, and it is shown how the effective diffusion coefficients of the sheath and the inter- 
stitial fluid can be obtained. 


Consider the situation in which a nerve trunk saturated with some sub- 
stance s is suddenly immersed in a solution free of this substance. We sup- 
pose that the amount remaining relative to the original amount after a 
time ¢ is a measured quantity G,;(#). It is the purpose here to derive an 
approximate expression for the quantity G(¢), which for the proper choice 
of the diffusion and permeability coefficients and the geometric and solu- 
bility factors will to a close degree of approximation equal the observed 
function Gz(t). 

Let D = Dgl where D is the coefficient and g and / are factors arising 
due to the reduced area for diffusion and increased path length (Schmidt, 
1952) due to the presence of the fibers. We shall assume that all fibers are 
equivalent with respect to the ratio of the permeability coefficient / to the 
fiber diameter p. Let the cross-section of the nerve, whose diameter is 2R, 
be composed of a fraction f, due to the fibers and a fraction fp = 1 — fa 
due to interstitial material. Similarly let a be the fraction of the actual 
volume of fibers which is available to the substance and let 6 be the frac- 
tion of the volume outside the fibers which is available to the substance. 

Case 1. Approximation method for the case of nerve without sheath for a 
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substance which does not penetrate the fibers appreciably or penetrates them 
very freely. Following the approximation method (Landahl, 1953), we sup- 
pose that the concentration EZ at some distance x from the center is equal 
to Eo, the initial value, for x < R—r, zero for x = R, and linear and continu- 
ous from R — r to R, the distance ry being a function of time. Since the 
integral of the flow rate is equal to the total amount which has left the 
nerve of length L, we may write 


= t 
2DaLEoR f O =  fybaLEor (3R—1). (1) 


The term on the left of (1) is the integral of the flow. The term on the right 
is the volume integral of Ey — £. 

Differentiating with respect to ¢ and setting D = D/df. (r = dr/dt), we 
find 


rr(3R—2r) =6RD. (2) 


On integrating we find, since r = 0 at ¢ = 0, 


r?(R—4r) =4RDt. (3) 
If we introduce 
Dt r 
ene ee Pee (4) 


then, since y < R, we may introduce for r in the parenthesis in (3) its first 
approximation. In expanding, we find an expression which is fairly accu- 
rate over the range y < 1 


y=2V7r+48r4+.... (5) 


The approximation (5) differs from (3) to the largest extent near r = .05, 
where y = .536 instead of .508. When 7 = 7, = 5/36 = .139, y = 1 or 
r = R in (3), while y = 0.99 from (5). 

After 4, = R’71/D, let E be approximated by a function which is 
linear with « being equal to fE) at « = 0 and zero at « = R, where f isa 
function of #. Equating the time integral of the flow with the total amount 
which has left the nerve, we find an expression for which the solution is 


f= e—8r—) 7 T>T71 = 5. (6) 


If Gi(r) = G(d), then to the same degree of approximation as above, 


eo 
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since G = 1 — y+ y?/3, 
Gi(r) =1-2V7r—.44474-237 24107, ron, ) 
7 
G; (7) = 4 e—6rt (5/6) 5 T= Tie f ; 


If when r(¢) becomes equal to some value R — RA, we suppose that r 
becomes constant but that the concentration within the range x < RA 
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FicureE 1. The fraction G; of the substance remaining in a nerve without a sheath plotted 
against r = Dt/R? for various values of 7’, the ratio of the initial amounts inside and outside 
of the fibers, and values of az = 2k2bf,R?/paglD, A being set equal to 3. The numbers, not 
otherwise specified, are values of a2. For the case in which y’ = az = 0, a comparison is given 
between the results (dotted) of the exact solution (Hill, 1928) and the approximation method 
(solid). 


decreases, being independent of x for x < RA, then one finds instead of (7) 


Gyula 1 — 3 7 — 4449 2.37 PY or <tT \ - 
Gi(r) =4 (1 FAA?) eVGA T) | or art, 


where 77 is given by (5) with y replaced by (1 — A). When A approaches 
zero, the above expressions reduce to (7). But even if A is one-half, the two 
expressions are almost the same for values of G down to a few per cent. In 
Figure 1, expression (8) with A = 1/2 is compared with the exact solution 


(Hill, loc. cit.). Expression (7) starts out with (8), approaching closer to 


the exact solution for longer times. The remainder of the curves shown in 
the figure will be discussed subsequently. The agreement among the three 
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expressions is seen to be satisfactory. However, it should be noted that the 
diffusion coefficient is somewhat overestimated by the approximation 
method if only the first part of the curve is used, whereas it is underesti- 
mated if only small values of G are used. 

If the permeability & for a particular substance being considered is 
large enough, e.g., k > 3pD/f.R?, then the above expression (7) or (8) 
gives the solution for this case also if bf, is replaced by bf, + af. wherever 
it occurs in those expressions. 

Case 2. Nerve without sheath. Permeability of fibers moderately small, 
k « 3pD/f,R?. If, in the process of equilibration, the total flow into the 
interstitial region is small compared with the amount that flows through 
that region, we may without appreciable error suppose that the fibers are 
confined to the central region of the nerve (cf. Shanes, 1951). Let AR be 
the radius of this region. In the region outside of A, we suppose that the 
effective volume occupied by the interstitial materials is the same as before. 
The volume occupied by fibers is treated as if unavailable for the sub- 
stance. Within A, we suppose that effective volume occupied by the inter- 
stitial materials is again unchanged but the volume occupied by fibers is 
that of the fibers in all the nerve. Hence for tr < r* the function G,(r) is 
given by (8) if one takes into account the amount within the fibers. We 
shall consider the case in which the permeability to flow into the fibers, 
ki, and the permeability to flow out from the fibers, ke, are not necessarily 
equal. Then k:Ey) = k2Co, where C» is the initial value of the concentration 
C within the fibers. Initially the ratio of the amount within the fibers to 
that in the interstitial fluid is given by 


ao afak, 
tame (9) 


The modified expression for G;(r) is given in terms of the parameter ’ by 
the first expression (14) below, for + < r*, the time at which 7 
= (1—A)R. 

After r has become equal to (1 — A)R, we proceed as follows. Let the 
concentration within the ring r = RA be treated as being independent of 
position but dependent upon time. Let E*(#) be the concentration in the 
interstitial space within the cylinder of radius RA. To simplify matters 
we shall suppose that the material which is outside of RA is effectively 
within RA, the effective volume being increased so as to make the con- 
centration equal to Eo at r = rj. The loss to E* is then due to diffusion 
proportional to DE*/(R — RA), and due to the inward permeability, 
proportional to k,£*; the gain is similarly proportional to A:C. Likewise 
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the net loss to C is proportional to k:C — k,E*. If we introduce 


6k; faR? 2k: b f,R? 
6;= < cet b \ 
pglD(i+A+a%)*’ “~~ pagiD’ eb0) 
we obtain 
= 
—(6+p2 5) Et +A, (11) 
dC 
7 PE ajH*—aC. (12) 
The ratio of the amount remaining to that initially present is 
a +35 f,(1+A+A*) E* 
G ae 3 ee 
ae SASS aL A pan ales ok 


Solving for C and E* in (11) and (12) for 8, = a considerably less than 
6 we find, on substituting into (13), that G,(z) is given by the second ex- 
pression of the following relation 


2M ot 4447 — 2.377 Nicolae 


Gi(7r) =1—-— ES, tin, 

14 
ie eee re At A%) ert An) > ot aa: 
al 3 +7) ees te 


If the a’s and ’s are not considerably smaller than 6, the equations be- 
come considerably more complex, the a’s and #’s both entering. Note that 
of these only az enters (14). In Figure 1 values of Gi(r) are graphed for 
various values of y’ and az. The largest value of az for each y’ is graphed 
from (14) even though (14) does not hold too accurately for such large 
values of ap. In making this estimate it is supposed that af./bf> is of the 
order of one. 

If one calculates the total flow out of the fibers per unit time at ¢ = ¢, 
and compares this with the rate of flow out of the nerve one finds, if 
A = 0, a ratio f.k:R?/pD which is equal to 3(h1/k2)(fa@/fob) (a2/6). The 
result does not vary appreciably with A unless A is greater than one- 
half. In those cases in which k; ~ ho, faa ~ fxb, the above ratio of flows is 
small if ap « 6 as was assumed to obtain (14). 

The parameters may be determined as follows for this case. Let an esti- 
mate of A be made, its exact value being unimportant. For short times 
(1 — G)? plotted against time gives D/(1 + 7’)? from the slope. From this 
value ¢; can be estimated. For long times, the last term in the second ex- 
pression of (14) is neglible (a2 « 6) unless ’ is small, in which case a 
knowledge of D from D/(1 + 7’)? is sufficient to account for the entire 
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curve except for the unimportant range G(t) < y’ <1. If 7’ is appreci- 
able, the plot of log G against ¢ — #; gives a linear relation with y'/(y’ +1) 
as the extrapolated value of G at ¢ = hy, the slope giving the value of he 
= h,/a, Rand p being known. But from 7’ and k2/a one obtains ki = faki/ bf. 
Thus the effective diffusion coefficient D and the effective permeability 
coefficients , and # can be determined. To obtain D, ki, and ke it is neces- 
sary to obtain estimates of a, bf>/fa, and gl/bf, independently. 

Case 3. Nerve with a sheath. Substance does not penetrate fibers or pene- 
trates very freely. Let the radius of the core be R and the thickness of the 
sheath eR. Let D and D, be the diffusion coefficients in the core and 
sheath respectively, each being corrected for increased resistance to flow 
due to various structures. Let the initial concentration in the sheath be 
Eo, that in the core being yo, y being a partition coefficient. Let v be the 
fraction of sheath volume available to the substance and let g = 2ve/y dfs 
be the fraction of the substance initially in the sheath. From the time 
zero until a time ¢,, we may calculate G(#) or G(r) from 


2 


q(Vr—.227+...) ig ole iets 3 


ante ae ere 4 


(15) 


where #, is given by (5) with y = e. In the above we have supposed e is 
not too large a fraction so that the true value e; should be replaced by 
€:(1 + de) to make the first approximation adjustment for the volume 
error resulting from dropping higher powers of e. 

To find a solution for ¢ > ¢,, we follow the same procedure as above. 
Equating the flow to the boundary at R to the flow away from the 
boundary, we find a relation which, solved for f, fEy being the concentra- 
tion in the sheath at R, may be written as 

Dey 
ed ie Ce 
For this purpose the concentration is taken to be yE» for x < R — # 
yfEo at « = R-, fEy at R*, and zero at « = R + eR. Since the time inte- 
gral of the flow equals the total amount which has left the nerve, we find, 
in a manner similar to that used to obtain (1), 


— rtfdt 
2eL EoD, f LA = wLER® (ve — fev-+ yb fey (3-=y) C1 19) MT) 


—— 
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Differentiating both sides of (17) and eliminating f by using (16), we find 
(y = r/R), if A = D,/Dye, 


[3g+12y +6 (A 1) y*—4Ay") 5 = 1214 Ay) B. (18) 
On integrating, (18) becomes 


2D (3A +1) G44) 
eat aa et aac, 2 tg BAD! 
(19) 
ee 
+2 mptar en, 
If y is obtained as a function of ¢ from (19), then G;(¢), with y as a parame- 
ter, is given by 
2 
Ay (1-y+%)4144¢ 
(1-9) Ay) ; 


Gili) = le<t<ke, (20) 


t. satisfying y(t.) = 1. 
After the time f2, when r = R, we approximate the concentration by a 
function made up of linear parts in which E(0) = FyEy, E(R-) = yfEo, 
E(R*) = fEo, and E(R + eR) = 0, f and F being functions of the time f. 
Using the same argument as above we find 
pat DeyF _ 

D,+ Dey 


C25) 
and 


- — wef + bvb fy EF2N IR. (22) 


Since F(t.) = 1, we find from (21) and (22) 


G() = 4S e—2AD(t—ty)/(1+1/2q+1/3 A) R? ; is ye AORN) 

The parameters may be determined from G(f) as follows. We suppose 
R and eR have been measured and that g, the ratio between the amount 
initially present in the sheath to that in the core, has been determined. 
For short times a plot of (1 — G)® against ¢ should be linear [from (15)] 
giving D,/v, denoted by D., from the slope. The quantity D, is then the 
effective diffusion coefficient of the substance in the sheath. For long times 
a plot of log G against ¢ should be linear according to (23), giving 
s = AD/(1 + 3g + 1/34) as the slope. Since we now know D,/v and 
also the slope s, then since s also equals (D,/v)q/2E(1 + 29 + 1/3A) we 
may find A from 


3(/D\ 49 
4=5(*) 4-3-4. (24) 
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With A now determined, the effective diffusion coefficient for the inter- 
stitial fluid, D, can be obtained from the expression for s. 

If expression (24) gives an ambiguous result for A, the following method 
may be useful. The intercept of log G versus ¢ — fz, extrapolated back to 
to, is (1 + 4g + 1/34)/(1 + A)(1 +) from which A can be calculated. 
However, fz is not known. Now #, can be obtained directly from (15) as 
the value of ¢ for which G(é) = G(é.) = (1 + 4¢)/(1 + 9). From (19), te 
satisfies the relation 


R75. 5 fp lenee vi 
sBiletis 4taatG a 54) log (1 + A) J. (25) 


ty =te+ 


Choose a value of & from which an estimated value A; of A can be ob- 
tained. Introducing this value into the expression for the slope gives D: 
an estimated value of D. Substituting these values of A; and D, for A 
and D in (25) should give the value of ¢, chosen. If the value obtained is 
not sufficiently close to that chosen, the process may be repeated for an- 
other value chosen to represent 2. If a satisfactory value is obtained in this 
way, the value of D is then known as well as the value of A which then 
determines D,/v or D,/ybf, through the definition of A. This should be 
compatible with the value obtained for short ¢’s. With D and A known, the 
branch of G(é) for t, < ¢ < f2 can be calculated readily through (19) and 
(20). Expression (19) now involves no parameters and can be plotted 
numerically. Similarly (20) involves no parameters. Agreement between 
theory and experiment in this range is meaningful since values in this 
range were not used to obtain the parameters. 

For the case in which the substance penetrates the fibers very freely 
we may simply replace bf, by fs + afay*, y* being a partition coefficient 
between the fibers and surrounding material. 

If the core is experimentally removed and replaced by a uniform liquid 
or gel, the above results hold, if the following changes and substitutions 
are made: bf = 1, D = D, where D, is the diffusion coefficient of the 
solute molecules in the introduced material, y is now the partition coeffi- 
cient between the substance of the sheath and the material replacing the 
core. 


The author is indebted to Drs. G. Schmidt and G. Karreman for read- 
ing and discussing the manuscript. 
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The attention of biologists is directed to the analogy between a metabolizing system and 
an electrical circuit with feedback and finite time delay. Periodicities can be explained on the 
basis of linear chemical kinetics and linear differential equations, although the latter are of 
infinite order. 


The purpose of this note is to suggest how periodicities in open systems 
may arise naturally from linear reaction kinetics. Previous explanations 
of periodicity have not been altogether satisfactory, as pointed out by 
F. E. Smith (1952) in his review of population theories, and more recently 
by J. Z. Hearon (1953). Even the classical Lotka-Volterra equations re- 
quire special assumptions. 

Another approach to the problem lies in striking an analogy between a 
biological system and a control system containing feedback. The compari- 
son is a direct one because metabolic reactions are known to involve not 
only chains but closed loops or cycles as well. N. Minorsky (1942) has 
shown how the presence of a finite time delay (metabolic lag) can cause 
steady oscillations where feedback exists. Consider a simple differential 
equation of the type which arises in a flow system 


gt+ax+ bi=0, (1) 


where x is any fluctuating property of the system and & signifies a retarded 
variable. Expansion of % in a Taylor series 


gr 
#=x(0—A0) =x—#0+%5—... + (— 1) (2) 


and insertion into (1) yields a linear differential equation of infinite order. 
Minorsky (loc. cit.) and T. Baron (1951) have shown that the solution of 
such a hysterodifferential equation has the form 

p= 4,62", (3) 
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where Z isa complex number and « is an initial disturbance to the system. 
The conditions for steady oscillation are 


ee ee TY 

Fa (4) 
and 

ea Tt 1=(5)- (5) 


The constant b in (1), (4), and (5) can be either positive or negative. 

The origin of metabolic lags has been discussed by C. N. Hinshelwood 
(1946) who points out that enzymatic reactions which are cyclically linked 
do not respond instantaneously to disturbances of the steady-state con- 
centrations. The arguments are straightforward and need not be repeated 
here. Interlocking chains and loops must be postulated in order to account 
for both time lag and feedback, but the reaction kinetics need not be 
complicated. 

An application of the above theory was reported recently by R. K. Finn 
and R. E. Wilson (1954). They found that the population dynamics of a 
yeast culture provided circumstantial evidence in support of the theory. 
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There are two kinds of scientific books worth reading. One is the monograph or treatise 
type, in which a more or less large field of science is presented in a systematic way, and in the 
form of a product as finished as possible at the given time. This kind of book may be consid- 
ered a source of knowledge then available. The other type of book may present a collection of 
chapters or individual articles which do not claim to be a complete and systematic treatment 
of the subject; however the reader not only finds interesting ideas there, but the reading as 
such suggests new ideas. Such books are useful. For, although a rough and unfinished idea per 
se does not even remotely have the value of a well-elaborated scientific study, yet no elaborate 
study, no important theory, can be developed without first having a few rough ideas. 

The book under consideration definitely belongs to the second category: it is a collection of 
essays. As the editor states in the Introduction (p. 2): “The papers in this volume are of a very 
different degree of maturity. They range from authoritative reviews of well-known facts to 
hesitant and tentative formulations of embryonic ideas.”’ He further states (p. 3): “We are 
aware of the fact that this volume is largely exploratory.” 

If the above is to be considered as a shortcoming, then the reviewer does not need to dwell 
on it, because the editor, and undoubtedly the authors, are fully aware of it, and duly warn 
the reader. If we evaluate the book from the point of view of how many ideas it suggests to the 
reader, then, at least so far as this reviewer is concerned, it must be considered a great success. 

The heterogeneity of the material, both in regard to the content and the length of the differ- 
ent contributions, makes a detailed paper-by-paper review somewhat difficult and hardly de- 
sirable. We shall therefore briefly mention only a few papers, without implying that the others are 
in our opinion in any way less interesting or less important. Different readers will undoubtedly 
be interested in different papers, according to their personal preferences. The reviewer felt par- 
ticularly interested in the contributions by Henry Quastler, ‘The Measure of Specificity’”’ 
(pp. 41-71); “The Specificity of Elementary Biological Functions” (pp. 170-88) ; the contribu- 
tion by Sidney M. Dancoff and Henry Quastler, “The Information Content and Error Rate in 
Living Things” (pp. 263-74); as well as the contribution of Herman R. Branson, “Informa- 
tion Theory and the Structure of Proteins” (pp. 84-104). 

In the first of the above contributions Quastler outlines a mathematical approach to the 
measurement of specificity of enzymes, by introducing specificity matrices. The rows of a 
specificity matrix correspond to different enzymes, the columns to different substrates. The 
elements 7;; are either unity, when the ith enzyme can react with the jth substrate, or are zero, 
if the corresponding enzyme and substrate do not react. Such matrices represent observable 
data. The author suggests expressions for the measure of specificity, keeping in mind that the 
conclusion as to specificity of an enzyme will always suffer some degree of uncertainty, due to 
the finite number of enzymes and substrates that can actually be observed. Possible generaliza- 
tions are suggested to other biological specificities, such as recognition. The second paper con- 
tains further similar considerations on specificity of biological functions, such as antigenic speci- 
ficities and allelic specificity. 

The third paper, by Dancoff and Quastler, discusses the fascinating but very difficult prob- 
lem of information content of an organism. Though the authors use very ingenious argu- 
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ments, they cannot arrive with any certainty even at the order of magnitude of the quantity 
sought. All they can do is to estimate the order of magnitude of the order of magnitude. They 
state (p. 270) that the information content of a man is not more than 10" bits and not less than 
105 bits. This does not look like much of a quantitative result. The reviewer, however, heartily 
subscribes to the authors’ statement that “this is an extremely coarse estimate, but is better 
than no estimate at all.”’ After all, Avogadro’s number was not at first computed with several 
decimals, as it is now known. The problem which the authors are facing is also a much more 
difficult one. 

Tt would seem that the larger figure may even be too low. Henry Linschitz, in a paper on 
“The Information Content of a Bacterial Cell”? (pp. 253-62), arrives at the figure 10" bits by 
a different method. It is to be regretted that this divergence between two papers in the same 
volume is not discussed in more detail. 

Dancoff and Quastler compare the obtained values for information content with the infor- 
mation content of an average printed page, 104 bits in their paper. On the other hand, Herman 
R. Branson, dealing with an entirely different subject, points out (pp. 84-85) that protein 
molecules are constructed very much as a message, “since they consist of some definite ar- 
rangement of about 20 amino acid residues.’’ All this suggests that estimates of information 
content of protein molecules as well as whole organisms may perhaps be made from the number 
of words necessary to completely describe in a given language the molecules or organisms. For 
some protein molecules such complete description may be in sight. For an organism this is 
much more doubtful. In that case it is perhaps better to reverse the argument. If we can com- 
pute the information content of an organism from biophysical considerations, then by com- 
paring it with the information available in language form (books or papers) about that par- 
ticular organism, we may estimate how much more information we must discover before we 
approach complete knowledge. . 

The book contains several interesting and highly informative papers, such as: ‘“The Protein 
Synthesis and Immunochemistry” by Felix Hourowitz (pp. 125-46), ‘““Genes and Antigens”’ by 
M. J. Irwin (pp. 147-69), and “The Control of Blood Sugar Level’? by Douglas A. Bragdon, 
Olga Nalbandov, and James W. Osborne (pp. 191-207), in which no quantitative discussion 
of information content is attempted. In fact, the word “information” is hardly mentioned at 
all. Any reader with a slight knowledge of information theory, however, will see what possible 
connections these papers have with the subject. Perhaps an editorial summary following each 
of these papers and pointing out the various information theory problems suggested in them 
would have been desirable, but we do not consider this an important point. 

A more serious thing, in the reviewer’s opinion, is the complete absence of contributions 
dealing with information theory and the central nervous system, which may be the field par 
excellence for the use of such a theory. Although no explicit reference to information theory is 
made in the well-known paper of W. McCulloch and W. Pitts (1943), the connection is quite 
obvious. This is made explicit in the systematic elaboration of the McCulloch-Pitts’ approach 
by J. von Neumann (1952). In his interesting book J. T. Culbertson (1950) discussed possible 
neural mechanisms for recognition of visual patterns, and particularly investigated the prob- 
lems of how greatly a pattern may be deformed without ceasing to be recognizable. The con- 
nection between this problem and the problem of distortion in the theory of information is ob- 
vious. The work of Anatol Rapoport and his associates on random nets, and especially on their 
applications to rumor spread (see the series of papers which appeared in this Journal during 
the past four years), is also closely connected with problems of information theory. The review- 
er (1950) has also made a small contribution to neuro-biological and bio-social aspects of the 
theory of communication. The paper by E. Reich (1951) also may be mentioned here. It seems 
to the reviewer that one or two papers dealing with the information theory aspects of the 
central nervous system would have added to the proper balance of the material in the book. 

To illustrate the statement made at the beginning of this review that this book suggests 
ideas, we shall end this review, perhaps somewhat facetiously, by a discussion of possible appli- 
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cations of information theory to writing books and book reviews. The purpose of a scientific 
book (we at least hope!) is to store and convey information in a given field. The purpose of a 
review is to convey information about a book. It is therefore legitimate to attempt a mathe- 
matical theory of writing books and to find the optimal conditions which make a book good. 
At first it may seem that the optimal conditions consist of maximizing the amount of informa- 
tion per page, that is, in minimizing the redundancy. But a certain amount of redundancy may 
not only be desirable, but necessary. When presenting a new subject to young students who 
have never heard of it, a judicious amount of repetition is good pedagogy. Giving an exact 
abstract definition and then illustrating it by an example already constitutes a logical redun- 
dancy. But how useful it frequently is! The minimum of redundancy that is found in some 
well-known and excellent mathematical books (nomina sunt odiosa!) occasionally makes those 
books difficult to read even for mathematicians. 

The optimum amount of redundancy is a function of the information and intelligence of the 
reader for whom the book is written. The analytical form of this function is to be determined 
by an appropriate mathematical theory of learning. Writing a book even in a field which be- 
longs entirely to the domains of Her Majesty the Queen of Sciences is, alas, still more an art 
than a science. Is it not possible, however, that in the future it may become an exact science ? 

Tf a reviewer’s information and intelligence are exactly equal to the value for which the 
book has been optimized, then he will perceive as defects in the book only deviations from the 
optimal conditions. His criticism will be objective and unbiased. If, however, the reviewer’s 
information and intelligence deviate in any direction from the value for which the book is in- 
tended, then he will perceive shortcomings which are not due to the deviation of the book from 
the optimum, but to the reviewer’s personal characteristics. He may also perceive some advan- 
tages in the same way. If in the society of the future every individual will be tagged, through 
appropriate tests, as to his information and intelligence at a given time, expressed in appropri- 
ate units, then a reviewer will be able to calculate the correction for his personal bias. These 
are fantastic dreams of today, which may become reality in the future. 

In the absence of the above-outlined theory this reviewer is not able to determine whether 
the following statement is due to his personal bias or not, but he feels strongly that the book 
under review should be highly recommended to all those interested in this field. 


N. RASHEVSKY 
Committee on Mathematical Biology 
UNIVERSITY OF CHICAGO 
March 16, 1954 
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